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Abstract: We study the chiral phase transition in a magnetic field at finite temperature 
and chemical potential within the Sakai-Sugimoto model, a holographic top-down approach 
to (large-A^c) QCD. We consider the limit of a small separation of the flavor D8-branes, 
which corresponds to a dual field theory comparable to a Nambu-Jona Lasinio (NJL) 
model. Mapping out the surface of the chiral phase transition in the parameter space of 
magnetic field strength, quark chemical potential, and temperature, we find that for small 
temperatures the addition of a magnetic field decreases the critical chemical potential 
for chiral symmetry restoration - in contrast to the case of vanishing chemical potential 
where, in accordance with the familiar phenomenon of magnetic catalysis, the magnetic 
field favors the chirally broken phase. This "inverse magnetic catalysis" (IMC) appears 
to be associated with a previously found magnetic phase transition within the chirally 
symmetric phase that shows an intriguing similarity to a transition into the lowest Landau 
level. We estimate IMC to persist up to 10^^ G at low temperatures. 



Keywords: 



Gauge-gravity correspondence, QCD, Chiral Lagrangians. 



Contents 



|T]. Introduction ^ 

|2|. General setup ^ 

|2.l| Context and brief summary of the model ^ 

2.2 Action on the flavor branes and equations of motion ^ 

|3|. Chirally broken phase ^ 

|3.1| Solution in the / ^ 1 approximation ^ 

3.2 Discontinuity in the density 12 



^ Chirally symmetric phase 14 



]J Zero-temperature limit and "Landau level" transition 15 



|5]. Chiral phase transition 19 



|5.l| Zero magnetic field 19 
5.2 Zero temperature ^ 



5.3 Zero chemical potential: when do chiral and deconfinement transitions split? 22 



5.4 General results: inverse magnetic catalysis and comparison to NJL ^ 
|6|. Conclusions |30 



|A| . Semi-analytic solution for / ~ 1 in the chirally broken phase 32 



p| . Semi-analytic solution for T = in the chirally symmetric phase 34 



[^ . Analytic approximation for zero-temperature "Landau level" transition 35 



|D| . Analytic approximation for zero-temperature chiral transition at small 



magnetic field 36 



1. Introduction 

Spontaneous breaking of chiral symmetry in a system of relativistic fermions is profoundly 
affected by an external magnetic field. A sufficiently strong, homogeneous magnetic field 
results in an effective dimensional reduction of the dynamics of the system. As a conse- 
quence, an instability with respect to condensation of fermion-antifermion pairs, i.e., with 
respect to the formation of a chiral condensate, occurs even at arbitrarily weak coupling. 
This is analogous to the Bardeen-Cooper-Schrieffer (BCS) mechanism of fermion-fermion 
pairing in a superconductor, where the effective dimensional reduction is achieved by the 
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presence of a Fermi surface. The enhancing effect of the magnetic field on chiral symme- 
try breaking has been termed magnetic catalysis (MC) and has originally been discussed 
in Gross-Neveu |, || and Nambu-Jona-Lasinio (NJL) g, |, |] models and in QED @. 
Recently, this effect has also been reproduced in holographic models with flavor branes 
subjected to magnetic flelds ||^, |8|, p^ ]. 

In this paper, we consider the chiral phase transition under the influence of a magnetic 
fleld B at finite temperature T and quark chemical potential //. Our main interest is QCD 
where chiral symmetry is spontaneously broken for sufficiently small temperatures and 
chemical potentials, and where the effects of strong magnetic fields may be observable in 
the chiral transition at small /x and large T (namely in ultrarelativistic heavy-ion collisions) 
and also at small T and large in astrophysical systems (compact stars). Indeed, it has 
been argued that extremely strong magnetic fields of up to ~ 10^® G occur in non-central 
heavy-ion collisions and up to ~ 10^^ G at the surface of magnetars (possibly 
even up to ~ 10^^ G in the interior ||l^). Our results may also be relevant for graphene 
1 14] which is under much better experimental control. Fermion excitations in graphene 
are effectively relativistic and MC manifests itself in a nonzero Dirac mass induced by 
electron-hole pairing [15, 16, 17], analogous to the constituent quark mass induced by 



quark-antiquark and quark-hole pairing in the QCD context. 

In both QCD and condensed matter contexts it is important to develop a strong- 
coupling description of MC. To this end we employ the AdS/CFT correspondence |l^, |l9|, 
|20| , 21], more precisely the Sakai-Sugimoto model ]^, |2^. This holographic model, based 
on type-IIA string theory, is, in a certain (albeit inaccessible) limit dual to large- A'^c QCD. 
In contrast to most other holographic models, it accounts for the full chiral symmetry group 
by realizing left- and right-handed massless fermions through Nf D8- and D8-branes in a 
background of Nc D4-branes. Moreover, the model has a confined and a deconfined phase, 
realized by two different background geometries. In the original version of the model, 
where the D8- and D8-branes are maximally separated in a compact extra dimension, the 
deconfinement and chiral phase transitions are identical and happen at a certain T for all 
values of fi and B (provided that any backreaction on the background is neglected). Here 
we are interested in a different limit of the model where the distance of the flavor branes is 
small and where a much richer phase structure is obtained. This limit can be understood 
as the NJL hmit of the model |2|, ||, ||]. 

The NJL model in its original form approximates the fermionic interaction by a point- 
like four-fermion interaction. It has been employed for the chiral phase transition in the 
presence of a background magnetic field at finite fj, and/or T in refs. ]27, 28, 30, 31, 
32, 33]. We shall find a phase diagram which shows a striking qualitative resemblance 
with some of the NJL results. In particular, we shall discuss that at finite chemical po- 
tential and not too large magnetic field the chirally broken phase becomes disfavored by 
increasing the magnetic field, in stark contrast to MC. We term this effect inverse mag- 
netic catalysis (IMC) and present a simple physical explanation, employing the analogy to 
superconductivity. Moreover, we shall discuss a discontinuity in the quark density ]34] for 
small temperatures in the chirally symmetric phase. Our result for this discontinuity, in 
particular the comparison of its location with respect to the chiral phase transition to recent 
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NJL results, supports its interpretation as a transition to the lowest Landau level. This is 



remarkable since with the exception of the Sakai-Sugimoto model |34, 35 1, Landau-level- 
like structures have been discussed in the AdS/CFT literature only in bottom-up scenarios 
1 36, 37, |3^, Here we shall see that the top-down approach of the Sakai-Sugimoto 

model suggests the presence of a lowest Landau level, but no further de Haas-van Alphen 
oscillations from higher Landau levels. 

Our calculation builds upon previous work within the Sakai-Sugimoto model in the 
presence of a magnetic field ||3^, 41, 42], and generalizes the results for the chiral phase 
transition in the T-/i plane at B = [43| and the T-B plane at /i = [44, ^] to the 



entire T-^-B space. (For recent discussions of the chiral phase transition in a magnetic 
field within other holographic models see for instance refs. 47, |4^, 49|.) 

The remainder of the paper is organized as follows. In sec. ^ we explain the geometry 
of our holographic setup, introduce our notation, and derive the on-shell action and the 
equations of motion. In sees. ^ and ^ we treat the chirally broken and symmetric phases 
separately and point out discontinuities in the quark density in both phases. The main part 
of the paper is sec. ^ where the chiral phase transition is discussed. After discussing several 
limit cases in sec. 5.1 - 15.3| (including a discussion of the split of chiral and deconfinement 



transitions in sec. |5.3D we present our main results in sec. ^.4| before giving our conclusions 
in sec. n. 



2. General setup 

2.1 Context and brief summary of the model 

We consider the Sakai-Sugimoto model for one fiavor, Nf = 1, in the deconfined phase. The 
corresponding background geometry is given by the ten-dimensional supergravity descrip- 
tion of Nc D4-branes in type-IIA string theory compactified on a supersymmetry breaking 



Kaluza-Klein circle [21|. Fundamental fiavor degrees of freedom are implemented by Nf D8- 
and D8-branes which are separated asymptotically by a given distance L in the compactified 



dimension |22, 23 1. Employing the probe brane approximation Nc ^ Nf, the background 
geometry will be fixed throughout the paper, while two qualitatively different embeddings 
of the fiavor branes account for the chirally broken and chirally symmetric phases. The 
U{Nf) gauge symmetries on the D8- and D8-branes are interpreted as left- and right-handed 
global symmetries of the dual field theory which lives at the 4-|-l-dimensional boundary 
of the ten-dimensional space (including the compact extra dimension, which needs to be 
small to arrive at an effectively 3-|-l-dimensional field theory). In the case of disconnected 
flavor branes the system is invariant under the full chiral group U{Nf)L x U{Nf)pi while 
connected branes in the bulk lead to the smaller symmetry group U{Nf)L+ji, see fig. |l]. 
This refiects the usual spontaneous chiral symmetry breaking of massless quarks. 

In the original applications of the model the separation of the fiavor branes L is max- 
imal, L = 'iT/M}<iyi, i.e., the D8- and D8-branes are put on opposite ends of the circle with 
radius M^^. Here, Mkk is the Kaluza-Klein mass which sets the mass scale below which 
adjoint scalars and fermions decouple from the dynamics of the dual field theory. For this 
maximal separation of flavor branes chiral symmetry is broken if and only if the system is 



-3- 



u 



D8 



D8 




U = QO 
U = Uo 



u = u 



Figure 1: Schematic picture of the chirally symmetric (left) and chirally broken (right) phases 
in the deconfined phase of the Sakai-Sugimoto model (in the confined phase, the X4-u subspace is 
cigar-shaped, not cylinder-shaped). The X4-direction is compactified on a circle with radius M^^, 
and u is the holographic coordinate with u = oo being the boundary where the dual field theory 
lives. If the distance L between the flavor branes is sufficiently small, a deconfined, chirally broken 
phase (right figure) becomes possible. In this phase, the connected flavor branes are embedded 
nontrivially in the background according to a function 14 (u) which has to be determined from the 
equations of motion which couple X4 to the gauge fields on the flavor branes. The location of the 
tip of the connected branes uq is part of this solution. In this paper, we work in the "NJL limit" 
of the model, where uq ^ ut, i.e., the tip of the branes is far away from the horizon. For a fixed 
temperature T - which fixes ut - the distance between uq and ut can be made arbitrarily large by 
choosing a sufficiently small separation L. 



confined. In other words, the chiral phase transition is dictated entirely by the background 
geometry. As a consequence, in the probe brane approximation the chiral transition is 
unaffected by all quantities that live on the flavor branes such as chemical potential and 
magnetic field. As an alternative to going beyond the probe brane approximation - which 
is very difficult -, such a "rigid" behavior can be softened by choosing a smaller separation 
of the flavor branes. This leads to a much richer phase structure, in particular a decoupling 
of the chiral and deconfinement phase transitions becomes possible. In fact, this decoupling 
is realized for values of L below a critical value L < 0.3tt/Mkk |5^) yielding a deconfined, 
chirally broken phase for sufficiently small chemical potentials [^]. Now that both con- 
nected and disconnected fiavor branes are possible solutions in the deconfined background 
geometry, also a magnetic field affects the chiral transition. It has been shown that for 
vanishing chemical potential, the critical temperature above which chiral symmetry is re- 
stored increases with increasing magnetic field in accordance with expectations 
from MC, as explained in the introduction. A simple consequence is that, in a certain 
regime of non- maximal separations L, a magnetic field may induce a splitting of chiral and 
deconfinement phase transitions. We discuss this effect in more detail in sec. |5.3| . Such a 
splitting has been observed in a linear sigma model coupled to quarks and Polyakov loop 
|51] and an NJL model with Polyakov loop (PNJL) (see however ref. |53|), but has not 
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been seen in lattice QCD calculations [^ ]. 

For a large separation of the flavor branes, i.e., L of the order of tt/MkK) the main 



features of the phase diagram in the T-^ plane are the same as for large- iVc QCD |55|. In 
the opposite limit L <^ iT/M}<i}^ the connected flavor branes are far away from the horizon so 
that the effect of confinement becomes less "visible" for the fundamental fermions. Hence 
we expect this limit to correspond to a field theory where the dynamics of chiral symmetry 
breaking decouples from that of the gluons. Therefore, by varying the distance L the Sakai- 
Sugimoto model interpolates between large- A'c QCD (L ~ vt/Mkk) and a (non-local) NJL 
model {L ^ tt/Mkk) where there are no gluons and no confinement [p4| , [25|] . Since in the 
former limit the many gluons dominate the phase diagram, the latter may in fact be an 
interesting limit for QCD at finite Nc, at least at comparatively low temperatures and high 
quark number densitites. Our main results correspond to the latter limit, which indeed 
will show many similarities to those obtained recently in NJL model calculations. 

2.2 Action on the flavor branes and equations of motion 

In the deconfined phase the induced (Euclidean) metric on the D8-branes is 



-j [f{U)dT^ + Sijdx'dx^ 



1 /[/^ ^ 



f{u) yR 



dU"^ + U^dnl J> , (2.1) 



where (r, xi, ^2, 2:3) are the coordinates of 3+1-dimensional space-time, X4 G [0,2tt/Mkk\ 
is the coordinate of the compactified extra dimension, R is the curvature radius of the 
background, and dQ^ is the metric of a four-sphere. The holographic coordinate on the 
flavor branes is denoted by U with U G [Ut, 00] (symmetric phase) and U S [[/q, 00] (broken 
phase), see fig. |l|, and 

f{U)^l-^, UT^i^Yr^R' (2.2) 



?73 ' ^ ~ V 3 7 

with the temperature T. We have used capital letters for the coordinates X4, U to reserve 
lower-case letters for their dimensionless versions introduced below. 

The action on the D8-branes has a Dirac-Born-Infeld (DBI) and Chern-Simons (CS) 
part, 

•S* = ^DBi + Scs ■ (2-3) 

Let us first discuss the DBI action. Its general form for one of the disconnected D8- and 
D8-branes in the chirally symmetric phase is 

d^x I dC/e-*Vdet(5 + 2TO'F) , (2.4) 
Jut 

where a' = with the string length is, the D8-brane tension Tg = (2Tr)~^£~^, the volume 
of the unit four-sphere V4 = 87r^/3, and the dilaton e* = QsiJJ / R)'^^^ with the string 
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coupling Qs- For brevity we have denoted the space-time integral by d^x although it is 
actually a Euclidean integral drd^x with imaginary time r, such that the integral over a 
space-time independent integrand (which is all we need in our calculation) yields V/T with 
the three-volume V. The DBI action for one half of the connected branes in the chirally 
broken phase is given by the same expression with the lower integration boundary Ut 
replaced by Uq. (In this general section, we shall give the expressions for the symmetric 
phase, but the broken phase is easily obtained via this simple replacement.) 

In our ansatz the only nonzero field strength components are F^o, Fus, B = Fi2- The 
field strength F12 is constant in the bulk and corresponds to a homogeneous magnetic field 
in the spatial 3-direction. Since the gauge symmetry on the flavor branes corresponds to 
a global symmetry at the boundary, B is not a dynamical magnetic field. However, this 
is not problematic in our context where we are interested in a fixed background magnetic 
field. 

We introduce the dimensionless quantities 

ao,3 = ~^^o,3 , b = 2TTa'B , = ^ ' ^ ^ ' 

where ^3 are the dimensionful gauge fields. Then, with the relation 

= TTgsNJl , (2.6) 
we can write the DBI action in the convenient form, also used in refs. 4C], 

5dbi = — [d^x r du yju^ + 62^,2^1 + _ ^,2 + ^3y^/2 _ (2.7) 

^ J J Ut 

Here the prime denotes the derivative with respect to ti, 

■^=6^(WF' ^^-^^ 

and 



nJ^ /" 47r \ ^ 

with the dimensionless temperature 



f{u) = l-^, ^T= f yj (2.9) 



t = TR. (2.10) 
The CS action is (in the gauge = 0) 

Scs = ^ I d^x I dU A^F^^Fp^e^'-'P'', (2.11) 



Ut 



where //, i^, . . . = 0,1,2,3 and £0123 = +1- Strictly speaking, this CS action is the action 
on the left-handed brane (D8-brane) . The corresponding action for the right-handed brane 
(D8-brane) has an overall minus sign. To avoid complications in the notation such as 
introducing left- and right handed gauge fields we shall only write expressions for the left- 
handed brane. This is sufficient for our purpose since we are mainly interested in the 
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free energy of the system, which does not distinguish between the left- and right-handed 
fermions. There are of course quantities, such as the currents, where the sign of the CS 
action becomes relevant. 

Within the above ansatz and using our dimensionless quantities, the CS part can be 
written as 

N f f°° r 1 

'S'cs = y J d'^x J du d2ai{aodua3 - a^duao) + ai{duaod2a3 - dua3d2ao) . (2.12) 

Here we have kept all terms that contribute to the equations of motion, although the two 
terms oc oi vanish in our on-shell action since 03 and oq do not depend on X2- 

If we worked with the action as given by eqs. and we would encounter an 



ambiguity in the currents: in the presence of a homogeneous magnetic field the currents 
defined via the usual AdS / CFT dictionary would deviate from the currents defined through 
thermodynamic relations. This problem was pointed out in refs. [40, In ref. [40| a 
modified action was suggested, 

5' = 5dbi + 5cs + A5, (2.13) 

with the additional contribution 

I (92 [oi {aoduas - as^^ao)] - du [ai (00^203 - 03^200)] | . (2.14) 



/\S = ^ [ (fx [ du- 



This term does not change the equations of motion since it is a boundary term, but it is 
not a usual holographic renormalization because the first term in the curly brackets is a 
term at the spatial boundary, not the holographic boundary. The new action S' is invariant 
under residual gauge transformations which do not vanish at the spatial boundary and 
removes the ambiguity in the currents. However, the correspondingly modified currents do 
not satisfy correct anomaly equations and reproduce the expected anomalous conductivities 



only up to a factor |42|. Here we are interested in the phase diagram, and not primarily 
in the anomalous conductivities, so we do not attempt to resolve this subtle issue; we 
simply follow the prescription with the modified action S' . The on-shell contribution of 
AS" is given solely by the term at the spatial boundary. This term becomes simply one 
half of the original CS part, and thus adding A5 effectively amounts to multiplying the 
original CS action by 3/2. Therefore, we expect our results to differ quantitatively, but not 
qualitatively, when we use the original action S instead. We can write our on-shell action 
for the left-handed flavor brane as 

VM 



™/2 
^4 



00 



+ ^^^&y duia^aQ- aoa'^) , (2.15) 
where we used d2ai = —b, and where the trivial space-time integral has been performed. 
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The equations of motion for oq, as, and X4 can be derived from eqs. ( ^I7| ) and ( ^.12 ), 



du 



du 



du 



IT 



360q 



3 ' 



72 



3ban 



; 



IT 



0. 



(2.16a) 



(2.16b) 



(2.16c) 



The left-hand (right-hand) sides of these equations originate from the DBI (CS) part of 
the action. Since the CS part does not depend on X4, the right-hand side of eq. ( |2.16c| ) 
vanishes. 

According to the AdS/CFT dictionary, the current is defined as 

dC 



J' 



dA' 



Consequently, the 0- and 3-components of the (left-handed) current are 



J' 



J' 



2-Ka' M 
2TTa' M 



-ba3 



2Tia' M 



-603(00) -|- c 



-6ao(oo) + D 



(2.17) 

(2.18a) 
(2.18b) 



^ 2 V ^ 2 7„^^ R 2 

where, in the second equality of each line, we have used the integrated form of the equations 
of motion ( 2.16 ) with the integration constants C and D. Due to our use of S", these 
results are identical to the ones which are obtained by taking the derivative of the free 
energy Q. with respect to the corresponding source. For example, is the charge density 
which is also obtained by the negative of the derivative of with respect to the chemical 
potential. Since VL will turn out to be very complicated in general, eq. ( 2.18a ) yields a 
simple alternative way to compute the density. 

In the subsequent sections we shall solve the equations of motion. For all /U, and b 
there are two classes of solutions. One with x'^ = 0, corresponding to straight, disconnected 
flavor branes and thus the chirally symmetric phase, and one with x'^ 7^ 0, corresponding to 
curved, connected flavor branes and thus the chirally broken phase, see fig. |l]. The general 
solution of the equations has to be found numerically, but we shall discuss various limits 
where semi-analytic solutions can be found. The solutions will then be inserted into the 
action in order to compare the free energies of the chirally broken and symmetric phase. 
This will lead us to our main result, the chiral phase transition as a critical surface in the 
t-^-b parameter space. 



3. Chirally broken phase 

3.1 Solution in the / — 1 approximation 

In general, the case of connected flavor branes is the more complicated one since besides 
the gauge flelds ag and 03 the equations of motion also contain the nontrivial function X4. 
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We simplify this case by approximating 



f{u) 1 (3.1) 

for all u on the flavor branes. This approximation is valid for sufficiently large uq ^> ut 
since for all u on the flavor branes we have u > uq, see fig. |l|. In principle, the approximation 
can - at a fixed temperature and thus fixed ut - be made arbitrarily good by decreasing 
the asymptotic distance of the flavor branes L. However, we have to keep in mind that we 
are interested in the critical temperature for the chiral phase transition. Suppose we choose 
L very small such that / ~ 1 is a good approximation at some small temperature. Then, 
increasing the temperature and keeping L fixed tends to invalidate our approximation 
because ut oc increases and approaches uq. But at some critical temperature the 
chirally symmetric phase takes over and thus our approximation only needs to be valid at 
temperatures below this (a priori unknown) critical temperature. At 6 = = 0, where the 
full treatment is simple, we have checked that our result for the critical temperature deviates 



by about 10% from the full result, see sec. 5.1. In the general case we have only solved the 



equations of motion in the limit / — 1, and thus have no quantitative comparison with the 
full result, but we have checked that the transition takes over before severe artifacts such 
as Uq < Ut occur in our approximation. Note also that within the approximation / — 1 
the broken phase becomes independent of T. (The chiral phase transition will still depend 
on T due to the T-dependence of the chirally restored phase.) 

With eq. (3.1) the integrated version of the equations of motion becomes 



a'Q^/u^ 


+ 62^2 






a'^Vu^ 


+ 52^2 


Vl + 42- 


Oq + u^x'l 


u^ x'^^Vu 


5 + 52^2 


Vi + 4^- 


Oq + u^x'l 



3603 + c , (3.2a) 
36ao + d, (3.2b) 



(3.2c) 



with integration constants c, d, and k. Our boundary conditions are 

a'o(uo) 



ao(oo) = /i, -f- — r = 0, (3.3a) 
03(^0) 

03(00)= J, «3('"o) = 0, (3.3b) 

x'iiuo) = 00 , 9 = / dux'4^, (3.3c) 

where we have introduced the dimensionless separation 

£ = I . (3.4) 

These boundary conditions arise as follows. First, we require the temporal component of the 
gauge field to approach the quark chemical potential // at the holographic boundary. Since 
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the chemical potential is the same for left- and right-handed quarks, uq must be symmetric, 
i.e., if ao is a smooth function along the entire connected branes, its derivative at the tip 
of the brane must vanish, aQ(uo) = 0. It turns out, however, that in the given choice of 
coordinates the more general version of this boundary condition is aQ{uQ) / a'^{uo) = 0. This 
takes into account that for finite magnetic field oq has a cusp at ti = uq, i.e., ao(n — )■ uq) 
approaches two values with the same magnitude but opposite sign depending on whether 
one approaches uq from the left or right. At the same time a'^iuo) becomes infinite. This 
apparent singularity can be removed by a coordinate change, for instance to the variable 
z used in ref. [|^], defined through u = (uq + uqz'^)^/^ (and vice versa, i.e., the smooth 



solutions in terms of z of ref. |42| acquire the same cusp in ao after changing the coordinate 
to u). 

For the spatial component in the direction of the magnetic field 03 we must allow for 
a nonzero value j at the holographic boundary which has to be determined dynamically 
by minimization of the free energy. It has been shown for the technically simpler cases of 



maximally separated branes and /or the Yang-Mills approximation of the DBI action [RS, 40 



41] that J assumes a nonzero value in the presence of a magnetic field. This corresponds 



to an anisotropic chiral condensate and thus, viewing this condensate as a superfluid, 
J corresponds to a supercurrent |41]. As a consequence, the system acquires nonzero 
baryon number, even for baryon chemical potentials smaller than the baryon mass 

Hill- 
Since the chiral condensate carries axial, not vector, charge, j is an axial supercurrent. 
Consequently, the boundary value at the other asymptotic end of the connected branes 
(where the right-handed fermions live) must be —j, leading to an antisymmetric gauge 
field 03 and thus to the boundary condition 03 (uq) = 0. 

Finally, the first boundary condition in eq. ( ^.3c| ) says that the connected branes "turn 
around" smoothly at u = uq while the second one says that the asymptotic (dimensionless) 
separation of the branes is i. 



We can solve eqs. (3.2) semi-analytically by generalizing the method introduced in 
ref. to the case x'^ 7^ 0. In this way, the differential equations can be reduced to two 
coupled algebraic equations which have to be solved numerically. We defer the details of 
this procedure to appendix The result can be written as follows. We introduce the 
constant rj via 

nf,= lim4M. (3.5) 



A nonzero r] implies that not only X4, but also the derivative of 03 becomes infinite at the 
tip of the branes uq. Moreover, we define the new variable y through 

, 3/2 J 

y{u)=3by^lT7f (3.6) 
Jun \/q{v) 



with 



g{u) ^ (77^ + l){u^ + _ ( ^2!^ ^ 1 I (^8 + ^2^5) _ (3 7) 
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In terms of these quantities, the solution for the gauge fields is 



ao(y) = IJ' + 



sinh yc 



-(coshy - coshyoo) 



sinh i/c 



■ sinh y . 



(3.8a) 
(3.8b) 



where ?/oo = y{u = oo), and the embedding of one half of the connected flavor branes is 
given by 

dv 



(3.9) 



The functions oq, as, X4, are written in terms of ^, b (which are the externally fixed physical 
parameters), j (which has to be determined from minimizing the free energy), and the 
constants uq, rj, which are functions of i, j, and b and given by the coupled equations 



i = uT\f4 + ^Q [ 

^ Ju. 



du 



3 



sinh yoo 36 ^/l + if 

The dependence on the separation i can be eliminated by rescaling 

b^i%, J^l^J, no ^ fuo . 



(3.10a) 
(3.10b) 



(3.11) 



(rj and y^o are invariant under rescaling with i.) Employing these rescalings and changing 
the integration variable u — >■ i'^u is equivalent to simply setting i = 1 in eq. ( 3.10a ). 

We could now proceed by solving eqs. ( p. 10 ) for all /i, 6, and j, insert the result into 
the solutions ( p.8| ) and (p.9|), these solutions into the on-shell action ( p. 15 ) and minimize 
the resulting free energy with respect to j. However, there is a simpler way to determine 
the supercurrent j. We recall that the total axial current J^^ is obtained by taking the 
derivative of the free energy with respect to the corresponding source. Here, j plays the 
role of that source and thus we conclude that j extremizes the free energy if J3 ^ 0, which 



implies, using eq. ( 2.18b ), 



J 



- tanh yoo ■ 



(3.12) 



This result can now be inserted into eq. ( 3.10b| ) which eliminates j from the numerical 
calculation. Written in this way, j is the same as in ref. |40|, but note that yoo is different in 
this reference. The reason is that there maximally separated flavor branes were considered, 
and thus the chirally broken phase was discussed in the confined geometry. 
The free energy is 

^ = 2 — 5'on_shell ) (3.13) 

where the factor 2 takes into account both halves of the connected branes and where 



'^on-shell given in eq. ( 2.15 ). After inserting the solutions of the equations of motion and 
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after some algebra we can write the free energy as 



w 



0,u = M I du , 



3^, / / 2 , n5+62^2 2 \ 



(3.14) 



tanh He 



where the result below the curly bracket eliminates j and has been obtained by using eqs. 
( 3.10b ) and ( 3.12 ). Inserting the rescaled quantities from eq. ( 3.11| ) into Qij shows that the 
free energy and the chemical potential scale as 



n 



(3.15) 



In the remainder of the paper we shall set ^ = 1 in all plots (except for fig. |7|) for convenience. 
The i dependence of all curves can easily be recovered with the rescalings (|3.1l| ) and (|3.15| ). 
Finally, we can compute the quark number density n. With eq. ( 2.18a| ) we find 



o 3 2TTa' , ^, 3 2'Ka' , ^, 
n = = -—^Mhj = -——Mbfitanhyo 



2 R 



4 R 



(3.16) 



(Recall that eqs. ( 2.18| ) only take into account one half of the branes, hence the factor 2 in 
the definition of the total density n.) 

3.2 Discontinuity in the density 



Although we can compute the supercurrent j directly from eq. ( |3.12[) , let us first discuss the 
form of the free energy as a function of j. Solving eqs. ( ^.lOD shows that there are parameter 
regions where there is a unique solution for the pair (uq , rj) and parameter regions where 
there are three solutions. This is reflected in the free energy shown in the left panel of fig. 
|2|. To obtain this plot we have renormalized Q by subtracting the vacuum contribution 



U" + 62^5 



+ 62 



(3.17) 



U'A 



where we have used that rj = for j = 0, see eq. (3.10b). The curves of the renormalized 
potential as a function of b show that there is a first-order phase transition where j is 
discontinuous and thus, due to eq. ( 3.16| ), also the baryon density n/Nc. The discontinuity 
is shown explicitly in the right panel where j is plotted as a function of b for three different 
values of fi. The full result has been obtained numerically, but we can easily find analytic 
approximations for small and large values for the magnetic field. For small magnetic fields 
and small j, eqs. ( |3.10D give 



Uo 



U0P2 



(3.18) 
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Figure 2: Left panel: free energy fiy as a function of the supercurrent j in the chirally broken 
phase for a fixed chemical potential fi and several values of the magnetic field h [in the dimensionless 



units given in cq. (2.5), and with i = 1 which is equivalent to using the rescaled quantities from eqs. 

and ( ini )]. Since j is dynamically determined from minimization of fiy, we see that for the 
given /i there is a first-order phase transition where j is discontinuous. Right panel: supercurrent j 
as a function of 6 for three values of /i. The solid lines are the full numerical result while the dashed 
lines are the linear approximations for small 6 according to eq. (3.20). The middle curve /i = 0.7 



corresponds to the potentials shown in the left panel. All three curves approach the asymptotic 
limit ] — /i/2 for large magnetic field, i.e., the upper two curves change their curvature to approach 
this limit on a 6 scale much larger than shown here. 
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Figure 3: First-order phase transition (solid line) in the fi-h plane for the chirally broken phase with 
6, /i in the same units as in fig. 0. The supercurrent (and thus the baryon density) is discontinuous 
across this line, as shown in fig. g. The dashed line is the chiral phase transition at T = (see sec. 
1^). On the right-hand side of the dashed line, the ground state turns out to be chirally symmetric 
(disconnected fiavor branes) such that the solid line is only relevant for a metastable state. 



with 



Pi 



i: 



du 



2^T 



. 16; 



. 16; 



(3.19) 
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We can thus approximate tanhyoo — Uoo and eq. ( 3.12 ) becomes 



16Pf 



fih (small 6). (3.20) 



This simple linear form for j is compared to the full result in the right panel of fig. For 
large magnetic fields, tanhyco — 1 and thus j approaches ^/2. 

In fig. m we show the discontinuity in the baryon density in the h-ji plane. As sug- 
gested from the right panel of fig. ||, the discontinuity is only present for sufficiently large 
chemical potentials. The first-order phase transition line terminates in a critical point and 
approaches the ^ axis for small magnetic fields. The figure also shows the chiral phase 
transition at T = 0, to be computed and discussed in sec. ^. On the right-hand side of this 
line, chiral symmetry is restored. Therefore, the discontinuity in the density only occurs in 
a metastable phase and is probably of little physical relevance. We shall thus not display 
it in the phase diagrams in the subsequent sections. 



4. Chirally symmetric phase 



The chirally symmetric phase has been considered in ref. |34] within the same setup as 
discussed here. Nevertheless we shall discuss some of the details of this phase before 
we come to the chiral phase transition. One reason is that we work at fixed chemical 
potential, while in ref. |34| the density was held fixed. Furthermore, we shall elaborate 
on a discontinuity in the charge density within this phase, which resembles a transition to 
the lowest Landau level. A physical understanding of this discontinuity will turn out to be 
useful in the comparison of our phase diagrams with NJL model calculations. 

In the case of disconnected flavor branes the integrated form of the equations of motion 



(|2.16D becomes 



a'n' 



yi + fa'i - a(f 



3ba3 + C, (4.1a) 
3bao + D, (4.1b) 



with integration constants C and D. Since the branes are straight, we have set x'^ = 0. 
In this sense, the equations are simpler than for the case of connected branes. However, 
now we cannot use the approximation / — 1 because the branes extend all the way down 
to u = ut, see fig. |l|. In this sense, the equations are more difficult than the ones for the 
connected branes. In general, we have to solve these equations numerically. Our boundary 
conditions are 

ao(oo)=/i, ao(nT) = 03(00) = 0. (4.2) 

As in the chirally broken phase, the value of oq at the holographic boundary is identified 
with the chemical potential. In contrast to the broken phase, the boundary value of 03 
vanishes because there are no Goldstone modes without spontaneous symmetry breaking, 
and thus there cannot be any supercurrent of these modes. We also require to vanish 
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at the horizon, which is a regularity constraint^ [^3|. For a^, there is a priori no condition 
at the horizon. Because of /{ut) = ao{uT) = 0, eq. ( [4.1b ) immediately yields 



D = 0. (4.3) 

(Provided that a'^{uT) is finite, which is true in all solutions we consider.) The numer- 
ical evaluation of eqs. ([4.1D can be done with the "shooting method": we consider the 
two differential equations as an initial value problem by imposing the initial values at the 
boundary n = oo according to eq. (^). Then we solve the equations by letting the gauge 
fields evolve from u = oo to u = ut for all C from an appropriately chosen interval. (It 
turned out to be useful to implement this procedure by promoting the ordinary differential 
equations to partial differential equations with the additional variable C.) Then we deter- 
mine the value(s) of C for which the gauge field oq is "shot" to its correct value at the 
horizon, ao(uT) = 0. (Since in some cases two of these values for C are very close to each 
other, it is more convenient to reparametrize C — )• 36/i coth z^o in the numerics, motivated 
by the zero-temperature solution, see below.) 

4.1 Zero-temperature limit and "Landau level" transition 

For T = 0, we have ut = and thus / = 1. In this case the equations ( [4.1D can be solved 
semi-analytically. The solution is (see appendix ^ for details) 

ao(z) = — ^ — sinhz, (4.4a) 
smh Zoo 

03(2:) = — ^ — (cosh 2: — cosh Zoo ) , (4.4b) 
smh Zoo 



with the new variable 



sinh 2:00 

where Zoo = z{u = 00) has to be determined numerically from the relation 



z(u) = 36 / "^"^ , (4.5) 

y5 + 52^2 + j3feM)^ 



.U> + 52^2 + 



sinh"^ 2, 

Inserting the solution ( |4.4| ) into the on-shell action ( 2.15| ) yields the free energy 



r°o n5 + 6V + i^^ 3 

Q|l(t = 0)=AA/ du , ^"'^ - -AA6/i^ coth Zoo. (4.7) 

,5 + 52^2 + _{|aOL 2 

sinh 2:00 







The quark number density is obtained from eq. ( 2.18a ), 



n = 2J^ = Z^^Mbfi coth Zoo . (4.8) 



^In Ref. jsj] it was argued that this regularity constraint needs to be abandoned in the case of an axial 
chemical potential. However, for an ordinary chemical potential as considered here, gauge invariance implies 
that this constraint is not a physical restriction. 



-15- 




Figure 4: Critical lines in the b-^ plane across which there is a discontinuity in the density of 
the chirally symmetric phase for three different (dimensionless) temperatures t (to recover the £ 
dependence, t has to be replaced by t£). For all nonzero temperatures the critical line ends at a 
critical point. We discuss in the text and with the help of figs. |^ and ^ that the critical line is 
reminiscent of a Landau level transition, i.e., a population solely in the lowest Landau level above 
the line and populated higher Landau levels below the line (no further transitions between the 
higher Landau levels are seen in our model). The thin (blue) dashed line is the same line as in fig. 
^ and indicates the chiral phase transition at t = (to be computed in sec. |^), i.e., on the left-hand 
side of this line the (blue) t — critical line has no physical meaning. The thick (black) dashed line 
is the analytic approximation ( [1.10 ) to the numerical t = result. 



One solution of eq. ( [4.q) is Zqo = oo- In this case, the density becomes 

n = ^Bfig , (4.9) 

where fiq is the dimensionful quark chemical potential, fiq = R/{27ra')fi, see eq. ( ^.5| ) for 
the corresponding relation for the gauge fields. 

The numerical calculation shows that in certain regions of the b, /x parameter space 
there are two additional nontrivial solutions for z^o- Also for nonzero temperatures, where 
we solve the differential equations purely numerically, one or three solutions are found. 

(1) (2) (3) (3) 

When we find three solutions < < z^J < Zod , where Zoo = oo for T = 0, the 
intermediate solution is never a global minimum of the free energy, while z2d and z!^ 
compete for the lowest free energy. Where the global minimum jumps from z2d to z!^ , a 
first-order critical surface appears in the t-fi-b parameter space. This surface is bounded 
by a critical line such that two-dimensional cuts through this parameter space, say at fixed 
temperature, show a critical line which, for nonzero temperatures, ends at a critical point. 
This is shown in fig. ^. For zero temperature, the critical line is given by the approximate 
critical magnetic field 

6^(t = 0) ~ 0.095 ^^/^ (4.10) 

This result is derived in appendix ^ and compared to the full solution in fig. ^. The ground 
state above this critical line is given by the solution Zoo = oo and thus the corresponding 
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Figure 5: Quark number density n as a function of the magnetic field b for several temperatures t 
at a fixed chemical potential /i — 0.5. There is a first-order phase transition at a critical magnetic 
field for sufficiently small temperatures. At t = 0, the density for magnetic fields above this critical 
value is exactly that of a non- interacting Fermi gas in a magnetic field, see eq. (|4.9|). 
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Figure 6: Number density for a system of non-interacting, massless fermions as a function of the 
magnetic field for several temperatures (all in appropriate units of the quark chemical potential 
/ig). The oscillations are caused by the successive population of Landau levels and are smeared out 
for large temperatures. This plot should be compared with the holographic result in fig. ||, which 
shows some features comparable to a Landau level structure. 



density by eq. (|4.9| ). Below the critical line the state with a nontrivial solution z^x, < oo 
(which depends on b and fi) has the lowest free energy. In this case, the density is more 
complicated. Only for b <^ bc{t = 0) , we find the approximate behavior n oc fj-^^"^, because 
in this limit z^o oc b/fi^^"^, see appendix |y. For nonzero temperatures, all solutions for z^o 
are finite and they continuously merge into each other for sufficiently small fi. 

In fig. ^ we show the density as a function of b for several temperatures at a fixed 
There are interesting parallels and differences to the case of free massless fermions in a 
magnetic field. The free energy of Nc non-interacting spin-^ fermion species of charge 1 in 



a homogeneous magnetic field B = is 

^^f.ee = -^f;(2-5.o)E^ r^ln[l + e-K^-^'^^)/^l , (4.11) 

where I labels the Landau levels. (In general, B has to be replaced by \qB\ in this expression, 
where q is the charge of the fermions.) The factor 2 — 6iq takes into account that the lowest 
Landau level (LLL) is occupied by a single spin degree of freedom, while all other Landau 
levels are degenerate with respect to both spin projections. The single-particle excitations 
are e^g,/ = Y^/cg + 2Bl^ where is the projection of the momentum on the direction of the 
magnetic field. The density follows immediately by taking the derivative with respect to 

NcB e- ^ /"°° f^^3 . / N /X 

nfree = V(2 - Oio) > s / — /^(efca,^ - s/^g) > (4-12) 

1=0 s=± -^^ 



where frix) = (e^/^ + 1) ^ is the Fermi distribution function. At T = 0, the distribution 
acquires a sharp Fermi surface and the density can be written as 



Nc „ , N,B 



nfree(T = 0) = ^S/ig + ^ _ 2BI . (4.13) 

1=1 



Here we have separated the contribution from the LLL which is populated for arbitrarily 
large B. The higher Landau levels / > are, for a given chemical potential, only populated 
for sufficiently small magnetic fields which is reflected in the upper limit of the sum over /. 

We plot nfi-ee as a function of the magnetic field in fig. ^. At T = 0, there are cusps 
in the density curve (i.e., discontinuities in the second derivative of the thermodynamic 
potential) which are caused by the Landau levels. Coming from large B, where only the 
LLL is occupied, contributions from higher Landau levels set in successively at each of 
these cusps. At small B, the sum over discrete levels can be approximated by an integral, 
and the result approaches the constant nfj-cc{T = B = 0) = NcHq/{37r'^) plus a highly 
oscillatory contribution with amplitude proportional to B^^"^. For arbitrarily small nonzero 
temperature the cusps are smeared out. The oscillatory behavior survives for small T and 
then completely disappears for large T. This is due to the smearing of the Fermi surface, 
i.e., at any nonzero T strictly speaking all Landau levels are occupied. 

We can summarize the comparison of our holographic result for the chirally symmetric 
phase to the particle picture as follows. 

• Zero temperature. - For large magnetic fields, the holographic density behaves exactly 
(i.e., all geometric constants of the model drop out) like that of a system of non- 
interacting fermions; this can be seen by comparing eqs. (|49|) and (p^) . In the 
particle picture, all fermions sit in the LLL in this limit. 

At a certain value of the magnetic field, namely B = //q/2, the non-interacting 
system starts to populate the first Landau level. This manifests itself in a cusp in the 
density curve corresponding to a second order transition, with infinitely many more 
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as B is lowered. In the holographic system there is instead a single first-order phase 
transition at the point where, coming from large B, the apparent LLL behavior 
ends. The critical value of B at which this transition happens cannot be directly 
compared to the one in the particle picture since it involves the geometric constants 
of the model such as the curvature radius R. In dimensionless quantities, this value 
IS i.e., it goes with a different power of /x than in the case of free 

fermions. In other words, the effective mass through the magnetic field seems to 
behave as B'^^^, not as B^/^. 

At small magnetic fields, the density in both systems becomes approximately constant 
in B, for free particles nfree oc /x^, while in the Sakai-Sugimoto model n oc /i^/^. 
Whereas the free fermion system shows an oscillatory behavior due to the Landau 
levels, the holographic result does not seem to know about Landau levels other than 
/ = 0. 

• Nonzero temperature.- While the cusps in the density of the ordinary fermionic sys- 
tem are smeared out at any nonzero temperature, the first order phase transition 
in the holographic result survives for small temperatures (the larger the chemical 
potential, the larger the temperature below which the discontinuity persists, see fig. 
^). Eventually, for sufficiently large temperatures, in both cases the density becomes 
monotonically increasing with increasing magnetic field, i.e., the transition in the 
holographic result disappears. 

5. Chiral phase transition 

Since the chiral phase transition has to be determined numerically in general, the next three 
subsections are devoted to some limit cases where the calculation is more transparent. 
These subsections also serve to discuss the / ~ 1 approximation in the chirally broken 
phase and the possible split of chiral and deconfinement phase transitions. 

5.1 Zero magnetic field 

In the chirally broken phase at vanishing magnetic field 6 = 0, the location of the tip of 
the connected flavor branes is given by eq. ( |3.18| ), 

IGvr 

uo{b = 0) = — 

We recall that here we have employed the approximation / — 1. In this case, we see that 
there is a unique solution for uq for any given £. This solution can become arbitrarily 
small. We need to ensure, however, that uq > ut in order to avoid the artifact of the 
flavor branes hanging farther down than they are allowed to by the geometry, see fig. 
|l[ With the result ( ^.1| ) and the definition of ut in eq. (^^) this condition is equivalent 
to t < ^^1/2 '^ {jq) /r {jq) — 0.173/^, which yields a temperature limit at 6 = for the 
applicability of our approximation. 



£(i) 



0.52 



(5.1) 
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In the full treatment, there is a critical value for the separation i above which there is 
no solution for uq (the branes must be disconnected then). For separations smaller than 



this maximal value there are in fact two solutions, one of which is unstable |58| and which 
approaches ut for ^ — t- (i.e., when the connected flavor branes are very close together 
they stretch down almost to the horizon). This unstable solution does not exist in the 
/ ~ 1 approximation, where the unique solution is an approximation to the stable solution 
of the full calculation. 

At zero magnetic field we have j = rj = and thus the free energy of the chirally 
broken phase ( ^.14|) becomes 

= 0)=Af r du u^/^ , ^ , (5.2) 



V 



with Uq given by eq. (5.1). 

In the chirally symmetric phase, the equations of motion in the b 
by setting b 



as = in eqs. ( |0 
when evaluated at u 



limit are obtained 
This yields a simple differential equation for ag, which. 



oo, relates the integration constant C to the chemical potential, 

C 



du ■ , 



(5.3) 



The free energy can be obtained from eq. ( p.l5| ). Using the equation of motion for oq we 
have 



Qn{b = 0)=J\f 



du = 



(5.4) 



The chiral phase transition is now obtained by finding the zero of the free energy difference 



An = Qn - 



(5.5) 



(While each of the free energies is divergent, their difference is finite.) Even in the case 
6 = 0, the zero of Ail has to be found numerically in general. Our result for zero (and 
nonzero) magnetic fields is shown in the next subsection in the lower panel of fig. |9|. For 
vanishing chemical potential we find the analytic result 



An{b = /i = 0) _ 2 
Jf ~ 7 



7/2 



which, using eq. (| 



yields the critical temperature 



t 



8/7 



0.14 



(5.6) 



(5.7) 



This critical temperature is close to, but still below the upper limit for our approximation 
discussed above. Our approximate value deviates from the full result by about 10% (see 
fig. 6 in ref. [^]). We can use our result to estimate for which separations L there is 
a deconfined, chirally broken phase. This phase occurs if Tc = tci/L is larger than the 
critical temperature for deconfinement Tc^deconf. = -^KK/(2vr). Consequently, the critical L 
below which a deconfined chirally broken phase exists, is Lc ~ 0.277r/MKK (compared to 
Lc ~ O.SIvt/Mkk in the full calculation [50|). 
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5.2 Zero temperature 

At zero temperature, we can compute the critical chemical potential for vanishing b as well 
as for asymptotically large b analytically. Since in our / — 1 approximation the chirally 
broken phase does not depend on temperature, the location of the tip of the connected 
branes uq and the free energy at 6 = t = are simply given by the results of the previous 
subsection, eqs. and ( |5.2| ). In the chirally symmetric phase, the value of the constant 
C at t = can be determined from eq. 



^2/5 



(5.8) 



The corresponding free energy is given by inserting this value and ut 
As a result, the difference in free energies becomes 



into eq. (| 



An{b = t = 0) _2 
J7 ~ 7 



)r(|) 



vr 



1/2 



(5.9) 



which yields the critical chemical potential 

, i6vrii/i4rr(3)r (6)15/7 
ficib = t = o) = ^ ^^^^ 



r(^) 



r 4 



16/7 



0.44 

1^ 



(5.10) 



Since at t = we have / = 1, this result is exactly the same as in fig. 6 of ref. p3(| . 

At asymptotically large magnetic field, sinhyoo diverges and thus eqs. ( 3.10bD implies 
rj = while from eq. ( p.l0a| ) we obtain 

2 

1.23 

(5.11) 



Uo[b oo) = — 2- 



£li) 



For the free energy in the chirally broken phase we insert yoo oo into eq. (|3l^ ), while 
in the chirally symmetric phase we use Zqo — t" oo in eq. ( [1.7[ ). Consequently, 



An{t = 0,6 oo) _ 6 
J7 ~ 2 



2^1/2 



r(f) 



9 



(5.12) 



which yields the critical chemical potential at asymptotically large b, 

5/2 

ndt = 0, 6 — )• oo) 



3^2 



r(f) 



0.43 



(5.13) 



Hence the critical chemical potentials at 6 = and 6—7-00 are almost identical. This is 
confirmed by the numerical result for arbitrary 6 which is presented in sec. 5.4 in fig. |^. 

It is important to specify that we compare the free energies of the two phases at a 
fixed value of the microscopic magnetic field B, not the externally applied field H. Had we 
fixed H, we would have had to perform a Legendre transformation of our free energy, as 
done in ref. |41]. In general, the physical context dictates which field must be held fixed. 
Here we are mostly interested in a comparison with previous NJL calculations, where B is 
fixed. 
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5.3 Zero chemical potential: 
split? 

In the chirally broken phase at /i 
T] = and 



when do chiral and deconfinement transitions 



0, eq. ( |3.12| ) imphes j = and thus eqs. ( 3.10 ) yield 



du ■ 



V 



(5.14) 



'uo ■a?l'^^/u^ + h^w" - (ng + 62^5) ' 

which is an equation for uq, to be solved numerically. One finds that uq increases mono- 
tonically with h and saturates at a finite value for asymptotically large b. This value can 
be computed analytically and is given by eq. ( ^.11 ). (Once we let 5 — cxd, the values of rj 
and uq at the minimum of the free energy become independent of //.) 

In eq. ( ^.Ij ) we have seen that the tip of the branes can be lifted by decreasing their 
asymptotic separation I. Now we see that a magnetic field has a similar effect: for a fixed 
separation I a magnetic field increases uq from the value ( |5.1| ) at 6 = to the value (5.11) 
for 6 = oo. Large values of uq tend to favor the chirally broken phase: see for instance 
eq. ([5. 71) which shows that decreasing the asymptotic separation i (and thus increasing 



Uq) increases the critical temperature tc- Therefore, a magnetic field seems to favor chiral 
symmetry breaking, which is in accordance with the expectation of MC. We will discuss 
this in more detail in the next subsection, where we show that one cannot naively transfer 
this expectation to the case of nonzero chemical potential. 

The free energy of the chirally broken phase ( 3.14 ) becomes with j = r] 







Ou(^ = 0) =M duu^/"^ 



+ 6V 



"0 



Y^u^ + h'^u^ - {ul + 62 



while for the symmetric phase we have 

O||(^ = 0) =M 



POO 

/ du \/u' 

J Uj' 



(5.15) 



(5.16) 



The phase transition for arbitrary b must be determined numerically. The resulting critical 
line for zero (and nonzero) chemical potential is presented in the next subsection in the 
middle panel of fig. |9|. Using the complete function /, this line has been computed in refs. 
1 44, 45 1 . Here we continue with an analytic result for asymptotically large b [for 6 = the 
result is given in eq. (|5.7])]. In this case, the difference of free energies becomes 



An{fi = 0,6 oo) _ 6 
If " 2 



.2 1/2 r (f ) 



2 

Ux 



(5.17) 



and thus, using the asymptotic expression for uq ( 5.11| ), the critical temperature is 

5/4 



tc{fj, = 0, 6 — )• oo) 



r(t) 



0.19 



(5.18) 



Again, we can now determine the critical value of the separation of the flavor branes below 
which the chiral and deconfinement phase transitions split. We find Lc ~ 0.387r/MKK- 
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Figure 7: Region in the b-L parameter space where the chiral and deconfinement phase transitions 
spUt (white area) and where they coincide (gray area) . For intermediate separations L of the flavor 
branes - between the two dashed hnes - a sufficiently large magnetic field b induces a split in the 
phase transitions. (While in all other plots it is understood implicitly that b is rescaled hy £ — L/R, 
we have written this rescaling explicitly here since L also appears on the horizontal axis.) 



This critical value is larger than without magnetic field. We can determine Lc for arbitrary 
b numerically to obtain a phase diagram with regions where the phase transitions coincide 
and where they don't, see fig. 0. 

The geometric picture is as follows. Suppose we set ^ = and choose the temperature 
slightly (infinitesimally) larger than the deconfinement phase transition (this transition 
is, in the probe brane approximation Nc ^» Nf, completely determined by Mkk and 
independent of all quantities on the flavor branes such as /i, b, and the separation i). It is 
now possible to choose the asymptotic separation of the flavor branes such that they can 
connect (meaning that the connected branes constitute the ground state of the system). 
To this end, the separation has to be sufficiently small. There is a regime of asymptotic 
separations (to the left of the left dashed line in fig. ^) where the flavor branes always 
connect, even for zero magnetic fields; there is another regime (to the right of the right 
dashed line) where the flavor branes never connect, even for asymptotically large magnetic 
fields; and there is a regime (between the two dashed lines) where the larger the magnetic 
field the farther apart we can put the connected flavor branes. 

Another way to read this figure is to consider the horizontal axis as a parameter that 
interpolates between different dual field theories. Large separations L correspond to large- 
Nc QCD, where the gluon dynamics becomes important. In this case, the magnetic field 
cannot induce a split of chiral and deconfinement phase transitions. Small separations L 
correspond to an NJL-like model, where the chiral and gluon dynamics decouple. While 
sufficiently small separations split the phase transitions for arbitrary magnetic field, there 
is an interesting intermediate regime where only a sufficiently large magnetic field induces 
a split. 

Although fig. 1^ has been obtained for vanishing chemical potential, the conclusions are 
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easy to generalize to all values of //. As we shall see in the next subsection, the critical 
temperature for the chiral phase transition is maximal for // = 0. Consequently, if and only 
if there is a split in the phase transitions at = there is also a split for a finite regime of 
nonzero /i. 

5.4 General results: inverse magnetic catalysis and comparison to NJL 

We can now compute the free energy difference between chirally broken and chirally sym- 
metric phases for all b, t, and The resulting chiral phase transition is presented in figs. 
|8| - |lO|, where figs. ^ and ^ are two-dimensional cuts through the three-dimensional phase 
diagram shown in fig. In fig. ^ and in the three-dimensional plot we show, in addition 
to the chiral phase transition, the "Landau level" transition discussed in sec. (In fig. 
m we have omitted this transition in order to keep the plots simple.) 

At the chiral phase transition line, the baryon number density increases from a purely 
topological contribution, which can be viewed as a stack of domain walls ||3^, to 
one that is carried by chirally symmetric quarks. For simplicity, our calculation does not 
take into account "normal" baryonic matter in the chirally broken phase, which in the 
Sakai-Sugimoto model can be represented by D4 branes wrapped on the 5*4 within the 
D8 branes. To get an idea about the possible onset of a normal baryon density we have 
computed the zero-temperature constituent quark mass which, in the setup with D4 branes, 
is b) = uo{iJ,, b)/3 |5^, |5^ (in the same dimensionless units as //), and have plotted 
the (thin dotted) line /i = no(/U,6)/3 in fig. ^. For fi > m{iJ,,b), an admixture of normal 
baryonic matter may occur in the broken phase. To obtain the actual transition to the 
phase where topological and normal baryonic matter coexist, and where it again ends, 
a consistent calculation including the effect of the baryon mass on the embedding of the 
flavor branes would have to be performed, which is however beyond the scope of this paper. 

The most interesting observations resulting from our present calculations are as follows. 

Inverse magnetic catalysis.- In fig. ^ we see that by increasing the magnetic field 
up to 5 < 0.2/ £^ at zero temperature and finite chemical potential, the chirally broken 
phase becomes less favorable. As discussed in the introduction, one might have expected 
a magnetic field to favor chiral symmetry breaking due to "magnetic catalysis" (MC). We 
term the observed opposite effect "inverse magnetic catalysis" (IMC). For larger magnetic 
fields, b > 0.2/ the ph ase transition line bends back and the magnetic field tends to 
favor the chirally broken phase, as expected from MC.^ Note that the two opposite effects 
occur on different scales of the magnetic field: the right panel of fig. § shows a large scale 
on which the phase transition line approaches its asymptotic value in accordance with MC; 
on this scale the opposite IMC at small magnetic fields is barely visible. 

The IMC becomes less pronounced for nonzero temperatures but exists up to t < 0.1, 
as we see in the upper panel of fig. ^. It manifests itself also in the middle and lower 
panels of this figure. For instance, in the middle panel we see a monotonically increasing 
critical temperature for /i = 0. For nonzero ^u, however, the critical temperature becomes 

^Had we restricted ourselves to isotropic configurations without a supercurrent, the phase transition line 
would show a somewhat enhanced IMC and then a weaker MC at large b, asymptoting to the smaller value 
M= ^Mc(i = 0,fe^oo). 
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Figure 8: Chiral phase transition (solid line) at zero temperature in the plane on a small 
(left) and a large (right) scale of the magnetic field. The first-order critical line divides the chirally 
broken phase (xSb, white area) from the chirally restored phase (gray area). It starts at the value 
( 5.1Cl|) for 6 = and approaches the value ( 5.13 ) for asymptotically large b (marked by an arrow in 



the right panel). The behavior in between is one of our main results since it shows that for finite 
chemical potential the presence of a small magnetic field disfavors chiral symmetry breaking. The 
dashed line marks the discontinuity of the quark density in the chirally symmetric phase. This 



transition is reminiscent of a Landau level transition, as indicated in the figure, see sec. 4.1 for a 
discussion of this line. The thin dotted line = mo(m, ^)/3 marks the potential onset of "normal", 
chirally broken baryonic matter to the right of this line (which is not included in our calculation). 



non-monotonic. There is even an intermediate range of /i, here shown for = 0.3, for 
which sufficiently cold matter is chirally broken at small and at large magnetic field, but 
not in between. 

Comparison with NJL calculations. - It is interesting to compare our results with cor- 
responding NJL calculations [27, 29, |3^, For instance, our fig. ^ and fig. 4 in ref. |29| 



show an amazing agreement in the chiral phase transition lines throughout the t-b-^ space; 
in particular, both results show IMC for moderate magnetic fields.^ In QCD, it is expected 
that part of the normal chirally symmetric phase is replaced by a color superconductor 



1 61]. Prom the results of the NJL calculation in Ref. one can read off that also in this 
case IMC is present for small temperatures. 

For t = 0, there are interesting similarities between our fig. ^ with fig. 4 in ref. 
as well as with fig. 2 in ref. |^^. Namely, in the left panel of fig. ^ we see that roughly 
at the point where the (dashed) critical line ends at the (solid) phase transition line, the 
latter strongly bends to the left. Such a structure is also seen in the NJL results, where 
the critical line marks the onset of the first Landau level. In the NJL model, more critical 
lines end at the chiral phase transition line, in principle one for each additionally occupied 
Landau level, giving rise to de Haas-van Alphen oscillations in the transition line. These 
additional lines and corresponding oscillations are absent in our approach, suggesting a 
separated "LLL" from a continuum of "higher LL's", as indicated in fig. |8|. 



claims that the Sakai-Sugimoto model does not agree with NJL results; but in this reference 
the CS contribution was ignored. 
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Figure 9: Transition between the chirally broken (xSb) and chirally symmetric (xS) phases in 
the b-fi, t-b, and planes for several fixed temperatures, chemical potentials, and magnetic fields, 
respectively (i.e., all panels are two-dimensional cuts through the three-dimensional phase diagram 
shown in fig. |T^. All lines are first-order phase transition lines. The quantities b, jj,, t are di- 
mensionless and related to the dimensionful counterparts by appropriate factors of (27r times) the 
string tension a' and the curvature radius R; moreover, we have set the asymptotic separation £ of 
the flavor branes to 1, the I dependence is recovered by replacing b b£^, fi -> t U. For 
simplicity we have omitted the "Landau level" transition lines shown in fig. @. 



Besides the absence of higher Landau levels, there are more differences to the NJL 
model. In particular, all our phase transitions are of first order. For the chiral phase 
transition this is obvious from the geometric point of view since there is a discontinuous 
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Figure 10: Chiral phase transition (large, blue surface) in the parameter space of magnetic field 
6, temperature t and chemical potential fj,. Two-dimensional cuts through this space (sometimes 
with larger b than shown here) are presented in fig. ^ Wc have indicated by "LLL" that, at t — 0, 
the quark number density in the chirally symmetric (xS) phase on the large-6 side of the (small, 
green) critical surface is identical to that of free fermions in the lowest Landau level. In general, 
this critical surface indicates a discontinuity in the quark number density. 



transition from connected to disconnected flavor branes. More physically speaking, there is 
a jump in the density across all phase transition lines shown in the plots. The NJL model, 
however, shows first-order and second-order phase transition lines in the chiral limit p9[ |. 

Discussion.- Why does dense (holographic) matter not behave as suggested by MC 
and rather shows a complicated mixture of MC and IMC? Naively, one might think that 
there are two "forces" acting on the chiral condensate: the magnetic field tends to effec- 
tively increase the particle-antiparticle coupling, as explained in the original works about 
MC m, ^ ^, ^, and thus work in favor of chiral symmetry breaking; the chemical potential 
tends to split particles from antiparticles which puts a stress on a particle-antiparticle pair, 
thus working against chiral symmetry breaking. And indeed, our results show instances 
where these two effects, considered separately, can be observed. Without ^u, the critical 
temperature is increased by B, see middle panel of fig. without B, the critical temper- 
ature is decreased by fi, see lower panel of fig. |9[ This suggests that if we increased B 
at fixed /i, only one of the "forces" would be at work, favoring chiral symmetry breaking. 
However, in fig. ^ we have seen that this is not true in general: for a certain range of 
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chemical potentials chiral symmetry is restored upon increasing B. 

To explain this effect, recall the analogy of magnetically-induced chiral symmetry 
breaking and superconductivity ||5| at weak coupling. In the case of a superconductor, 
conventional BCS Cooper pairing between (massless) fermion species whose Fermi surfaces 
are split by a mismatch Sfj, is possible if the pairing gap A is sufficiently large. One can 
picture this situation as follows. Start from two different filled Fermi spheres whose radii 
differ by 5fi. For conventional Cooper pairing to happen at zero temperature, the Fermi 
surfaces must coincide. To this end, force both Fermi surfaces to the common, average 
Fermi surface /i. Creating this fictitious, intermediate state results in a free energy cost 
oc ji^SiJ? . But now pairing yields a gain in free energy oc A^^^. Consequently, pairing is 
possible if A is large compared to 5^ and breaks down otherwise. Working out the correct 
prefactors, one finds that Cooper pairing breaks down for Sj^L > A/-v/2, which is called the 
Clogston-Chandrasekhar relation |62, 63 1. 

We can transfer this picture to the chiral condensate in a strong magnetic field as 
follows. First we note that we can restrict ourselves to the physics of the LLL, since both 
NJL and holographic results show the strongest IMC in a regime where the higher Landau 
levels are empty. While for the usual superconductor at (5/i = the fermions that "want" 
to pair sit on the two-dimensional surface k = kp oi the Fermi sphere {kp being the Fermi 
momentum), the LLL fermions and antifermions that "want" to form a chiral condensate 
both sit, at // = 0, on the two-dimensional plane k^ = perpendicular to the magnetic field 
in the 3-direction. Now we switch on /i, which in our context is the analogue of (^/i because 
it separates the fermion surface from the antifermion surface. As above, we imagine to 
force the two planes back to their /i = position. The resulting energy cost is oc Bfi^. 
This can be seen from the zero-temperature limit of eq. ( 4.11 ) which shows that the LLL 
contribution is = —NcBfj?/{47r'^) which becomes $7 = for /i = 0. Independent of the 
precise form of the energy gain due to the formation of a chiral condensate - which is also 
expected to increase with B - our first important observation is that the cost increases with 
B. (This is almost as if, in the case of the superconductor, both and A increase upon 
increasing a single parameter.) Therefore, the competition between the effects of /i and 
B is more complicated than naively expected, and we understand why IMC can happen. 
Whether it does happen depends on the coupling strength as we now explain. 

In the weak-coupling limit, the free energy difference obtained in an NJL model is 
|6|,||, 

An(xB I — - fi^\ , (5.19) 

where M is the -B-dependent constituent quark mass. In this case, there is an exact analogy 
to the Clogston-Chandrasekhar relation, namely > M/-v/2, and increasing B at fixed fj, 
can only increase, never decrease, A^l (since M increases with B at weak coupling due to 
MC). This shows that IMC is not possible in the weak-coupling limit. Interestingly, the free 
energy difference in our holographic calculation assumes the same form for asymptotically 
large B if we identify M = uqR/ [lira') |5^, as we can see from eq. ( 5.12| ). Therefore, 
we might speculate that the limit of asymptotically large magnetic fields, where we do not 
observe IMC, is in some sense equivalent to the weak-coupling limit. The reason might be 
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that the magnitude of the constituent quark mass can be interpreted as a measure for the 
couphng strength, and in the given hmit the constituent mass (squared) is much smaller 
than B. However, as eq. ( ^.12| ) shows, the relation between the condensation energy and 
the constituent quark mass involves a more complicated numerical factor compared to the 
NJL model. (For a comparison with eq. ( 5.19| ) we need to consider an isotropic condensate, 
i.e., switch off the supercurrent; even after this modification the prefactor is different.) 

Our observation of IMC at smaller magnetic fields suggests that the free energy dif- 
ference must change qualitatively. The simplest way to see this is to use the small- -B ap- 
proximation for the chirally broken phase and the "LLL" result for the symmetric phase. 
From eqs. (13) and (|]2D we find with uq (x M 

An ~ const X - fi^B , (5.20) 

i.e., the condensation energy has dramatically changed while the cost of forming a conden- 
sate has remained the same. For small magnetic fields we have M = const + 0{B^) (see 
appendix ^ for the precise form of this expansion), and thus at some value of B the cost 
exceeds the gain, resulting in IMC. Two comments about eq. (5.20) are in order. Firstly, 
one might question the validity of this free energy because we have used an expansion in 
small B, although the "LLL" phase exists only at sufficiently large B. Indeed, to obtain 
a good approximation to the full result, at least 0{B'^) terms have to be included, see ap- 
pendix However, this does not change the conclusion regarding IMC whose qualitative 
form is well captured by the above linear approximation. Secondly, eq. ( ^.20| ) shows that 
our conclusions for a truly four-dimensional field theory have to be taken with some care: 
M''/^ does not have mass dimensions of a free energy density; the constant contains the 

1/2 

dimensionful factor Mj^j^, reflecting the extra dimension in our model. Nevertheless, the 
qualitative agreement with NJL calculations suggests that our observation is of general 
nature and may thus also be relevant for QCD. 

What is the range of magnetic fields in physical units for which IMC occurs? For a 
rough estimate let us match the 6 = values of our critical temperature at /i = and our 
critical chemical potential at t = to the approximate values from QCD, Tc — 150 MeV 
and (iq^c ~ 400 MeV (the former is in fact a cross-over rather than a critical temperature, as 
known from lattice calculations |6^, 37 1, while for the latter we only have a comparatively 
rough idea, see e.g. refs. |68, ^). We can express the dimensionful quantities as 

where we have reinstated the electric charge q. This shows that, expressed in terms of the 
dimensionful model parameters E? /{2TTa') and L, the different scalings with respect to I 
used in our plots arise naturally.^ The scale for the magnetic field is now found with the 
help of these relations, 

l*^^' 5.1 X 1019 G f ( , (5.22) 



bi^ ' V 400 MeV/ V 150 MeV 



*With ~ jv^MK'it recover the parameters k and Mkk, whose values are matched to the physical 
pion decay constant and rho meson mass in ref. note, however, that we cannot simply use these 

numerical values since they are only meaningful for a maximal separation L — — . 
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where we have inserted our results for tcf- and Hc^-"^ from eqs. (|5.7|) and ( ^.lOP , respectively. 
Now we can read off from fig. ^ that, at zero temperature, IMC occurs for magnetic fields 
up to \qB\ < 1.0 X 10^^ G and leads to a reduction of the critical chemical potential from 
(the matched value) ~ 400 MeV down to ~ 230 MeV. The phase that we have identified 
with the transition into the lowest Landau level occurs for \qB\ > 1.0 x 10^^ G. 

Finally, let us elaborate on the comparison to the NJL model calculations. We have 
seen interesting similarities in the results. However, in the NJL works mentioned so far 



27, 29, Bfl] only isotropic chiral condensates have been considered. As a consequence, the 



chirally broken phase has, at least at small chemical potential, vanishing baryon number. 
Only if the quark chemical potential becomes larger than the constituent quark mass, the 
baryon number may be nonzero in the broken phase [27, 30 1. This is different in our holo- 
graphic calculation. Here we have an anisotropic chiral condensate throughout the chirally 
broken phase, which manifests itself in the nonzero supercurrent j and a nonzero topo- 
logical baryon number. As discussed above, we have not included "normal" homogeneous 
baryonic matter into our calculation. In view of these different kinds of baryonic densities 
it appears quite remarkable that the holographic and NJL phase diagrams look similar. 

In one recent NJL calculation [31], however, a more general ansatz has been considered. 
And indeed, as in our calculation, an anisotropic chiral condensate is found to be favored 
throughout the chirally broken phase. Curiously, the resulting phase diagram in fig. 1 of 



ref. looks less similar to our result, compared to the phase diagrams in refs. |27, pO| , 
which are obtained with an isotropic condensate. We can only speculate whether this might 
be due to the specific choice of the coupling constant. For a more reliable comparison to the 
NJL phase diagram it is crucial to extend our holographic results by including "normal" 
baryonic matter. 

A further possible complication is the so-called "chiral shift" in the symmetric phase 
1 64, 7C]. NJL model calculations predict this difference in the dispersions of left- and 
right-handed fermions in the presence of a magnetic field, and it remains to be seen how it 
influences the chiral phase transition. In our holographic calculation the chiral shift would 
correspond to a nonvanishing boundary value of 03, which is absent in the symmetric phase. 



6. Conclusions 



We have discussed the effect of a magnetic field on the chiral phase transition in the de- 
confined phase of the Sakai-Sugimoto model. In the probe brane approximation Nc ^ Nf 
applied here, this chiral transition exists only under certain conditions. It does not exist 
if the asymptotic separation of the flavor D8- and D8-branes in the compactifled extra 
dimension is sufficiently large because then chiral and deconfinement phase transitions 
coincide, independent of the magnitude of the magnetic field. We have identified an in- 
termediate region for the separation where a magnetic field induces a split of the two 
transitions, allowing for a chirally broken, deconfined phase. For even smaller separations, 
chiral and deconfinement phase transition are distinct even for vanishing magnetic field (at 
zero chemical potential). This is the regime of separations we have considered in the main 
part of the paper. It corresponds to an NJL-like model on the field theory side since the 
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chiral dynamics completely decouples from the gluon dynamics, and confinement becomes 
irrelevant. In this sense, our results are not of direct relevance to QCD, at least not for 
the interplay of chiral symmetry breaking and confinement. They may still be relevant for 
qualitative features of the chiral phase transition in QCD, in particular at large chemical 
potential and low temperature, where gluonic degrees of freedom are less dominant at fi- 
nite Nc- Moreover our results may be of general interest for other (effectively) relativistic 
systems with flavor symmetry breaking in a magnetic field, for instance graphene. 

We have computed the critical surface of the chiral transition in the three-dimensional 
parameter space of temperature, chemical potential and magnetic field. The most interest- 
ing result is observed for small temperatures and a certain intermediate range of chemical 
potentials. In this case, starting from a chirally broken phase at zero magnetic field, a 
small magnetic field induces symmetry restoration before chiral symmetry is broken again 
at large magnetic fields. The tendency of a magnetic field to favor chiral symmetry breaking 
is well known from the so-called magnetic catalysis. Our observation for dense holographic 
matter is more complicated because for small magnetic fields we see the opposite effect, 
which we have termed inverse magnetic catalysis. We have explained this effect in a simple 
analogy with a superconductor with mismatched Fermi momenta. The essence of this ar- 
gument is that the magnetic field not only enhances the fermion-antifermion coupling and 
thus the energy gain from forming a chiral condensate but also enhances the free energy 
cost needed to form antifermion-fermion pairs in the presence of a chemical potential. 

We have also pointed out parallels and differences of our results to previous NJL 
model calculations. Most nontrivial features of our results, such as the inverse magnetic 
catalysis in certain regions of the phase diagram, can be observed in an NJL model as 
well, supporting the interpretation of this specific limit of the Sakai-Sugimoto model as 
a holographic strong-coupling, non-local version of NJL. In accordance with Ref. we 
have found that the Sakai-Sugimoto model shows indications of a Landau level structure. 
This is suggested by the dependence of the quark density on the magnetic field which we 
have compared in detail with the corresponding density in a usual particle picture. There is 
a first-order phase transition within the chirally restored phase whose location with respect 
to the chiral phase transition is comparable to that of the (second-order) transition into the 
lowest Landau level in the NJL model. However, in the Sakai-Sugimoto model there are 
no further de Haas-van Alphen oscillations and no additional transitions corresponding to 
the higher Landau levels of the particle picture, which seem to be replaced by a continuum 
of states in the holographic model. 

Our work opens several directions for future projects. A straightforward extension 
is to take into account the full effect of the curved geometry which becomes important 
for large temperatures (we have considered the f{u) ~ 1 approximation in the chirally 
broken phase which is a considerable technical simplification). Moreover, one should include 
"conventional" baryonic matter in the chirally broken phase, in addition to the baryon 
number induced by an anisotropic chiral condensate. We have discussed that such an 
extension is interesting in view of a more detailed comparison to the NJL model. It would 
also be interesting to get a deeper understanding of the apparent Landau level structure 
in the Sakai-Sugimoto model, in particular it would be important to understand in which 
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sense this structure is a strong-coupling version of the usual discrete Landau levels. 
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A. Semi-analytic solution for / ~ 1 in the chirally broken phase 



Here we explain the solution of the equations (^]^). We first take the pairwise ratio of 
these equations, 



a'o 8603 + c 

a'g 36ao + d 

q'q Sbas + c 
u^x'^ k 

a'g 3bao + d 
u^x'a k 



(A.la) 
(A.lb) 
(A.lc) 



Evaluating eq. ( [A.la ) at the point u = uq yields [using the boundary conditions (|3.3|)] 

c = 0. (A.2) 
With the definition of r] in eq. ( |3.5D , eq. ( 3.2c ) at n = -uq yields 



k 



(A.3) 



Here we have used that, if 77 7^ 0, we must have a'^{uo) = 00 due to the boundary condition 
2^4(^0) = 00 and the definition (3^). Now we rewrite eq. ( [A.la ) as 35aoag + dog = 360303 
and integrate this equation to obtain 



^bal + dao = ^baj + k, 



(A.4) 



with another integration constant k. Evaluating this at ti = 00 and using the boundary 
conditions in eq. ([3.3|), yields function of d, 



3 



(A.5) 



We shall return to the determination of the integration constants below. Now we first 
rewrite the differential equations. To this end we insert u^x'^ from eq. (A.lb) into eqs. 
( 3.2a ) and ( 3.2b[) . The resulting equations are equivalent to 



dyttQ — 03 , 
dya^ = ao + 



36 



(A.6a) 
(A.6b) 
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and a condition for the new variable y, 



3b 



(A.7) 



(Remember that the prime always denotes the derivative with respect to u] derivatives 
with respect to y are written explicitly.) We shall discuss the solution of eqs. (A. 6) below; 
first we use the result for y' to obtain the solution for X4 and more relations between the 
integration constants. From eq. ( A.lb| ) we obtain 



Consequently, the condition that X4 diverge aX u = uq implies that y' also has to diverge 
sX, u = uq. Hence, the denominator on the right-hand side of eq. ( [A.7] ) has to vanish for 
u = uq which implies 



1 + 



and thus y{u) is given by 



with g{v) defined in eq. (| 
solution for X4, 



J u 



"0 yg{v) 



(A.9) 



(A.IO) 



From eqs. (A. 8) and ( A. 3D we immediately obtain the final 



X4(n) = v?J'^^Jul + h^ul 



dv 



.3/2^ 



(A.ll) 



Inserting this expression into the boundary condition for the asymptotic separation i, see 
eq. (^.3c|) , yields one of the two equations relating rj and uq. 



i 



3/2 



du 



The differential equations (A. 6) are solved by 



d 



(A.12) 



ao(y) = ci coshy + C2 sinhy - — , 

36 

03(2/) = ci sinhy + C2 coshy . 
The boundary condition 03(^0) = becomes 03(2/ = 0) = which implies 

C2 = 0. 



(A.13a) 
(A.13b) 



(A.14) 



Then, with ?/oo = v{u = 00) we have the boundary conditions ao(2/oo) = and 03(2/00) = J 
and thus 

, rl 

3 



ci 



cosh yoo sinh yc 



(A.15) 
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Consequently, we arrive at the final solution for the gauge fields, 

3 



ao{y) 



sinh Uc 



-(coshy — cosh t/q 



sinhyc 



sinh y . 



(A.16a) 
(A.16b) 



The second remaining equation to determine uq and t] is obtained by rewriting eq. ( |3.5| ) 
with the help of eq. (|A.8|) , 



3bu', 



3/2 



V 



3bu, 



3/2 



y=o 



k sinh 



where in the last step we have used eq. ( 3.8b ). Inserting k from eq. (| 



yields 



sinh Uc 



3b 



(A.17) 



(A.18) 



As a check, this relation can also be obtained as follows: one solves eqs. (A. 5) and ( [A.9[ ) 
for d and inserts the result into eq. ( A.15| ). After a little algebra one arrives at eq. ( |A.l^ ). 

The complete solution of the equations of motion is thus given by eqs. ( A.16| ), ( |A.11| ), 
( A.12|), and (A.18|) which are all given in sec. 0| of the main part of the paper. 



B. Semi-analytic solution for T = in the chirally symmetric phase 



In this appendix we derive the solution (4^) to the equations of motion in the chirally 
symmetric phase at T = 0. 

Dividing eq. ( [4. la)) by eq. ( [4.1b| ) yields 

«o 36a3 + C 

~T = — ' i^-l 

ffg 6bao 

where eq. ([4.3|) has been used. We write this as SbagaQ = 3ba'.^a^ + Ca^ and integrate to 
obtain 

hal = hal + Cas + K, (B.2) 

with an integration constant K. This constant is easily determined by evaluating eq. (|B 
at u = oo, 

K = hf^\ (B.3) 



With the new variable 



z{u) = 3b / , , (B.4) 

Jo ^V^ + 62i;2 _ (35^)2 + ^2 



we can write the equations of motion (^^) as 



(J 

dzao = "3 + ^ , (B.5a) 
dzas = ao , (B.5b) 
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which are solved by 

ao{z) = Ci cosh z + C2 sinh z , (B.6a) 

C 

03(2) = Ci sinh z + C2 cosh z -. (B.6b) 

ou 

The integration constants Ci, C2, C can now be determined from the boundary conditions 
( |4.2| ), and the resulting solutions become 

oo(z) = — 7- — sinhz, (B.7a) 
smh Zoo 

osfz) = — ^ — (cosh z — cosh Zoo) , (B.7b) 
smh Zoo 



with Zoo = z{u = 00), and 



dv 

z{u) = 3b ==. (B. 

Jo + 52^2 + J^aOL 

V smh Zoo 



From the equation for z(n) WG obts-in an implicit Gqucition for ^oo? 



Zoo = 36/ , (B.9) 



sinh 2o 



The complete solution is thus given by Eqs. (B.7) with z(u) given in Eq. ( |B.8| ) and Zoo to 
be determined numerically from Eq. (|B.9|). 



C. Analytic approximation for zero-temperature "Landau level" transi- 
tion 



Here we derive eq. ( 4.1C| ), which is an approximation to the full numerical result for the 



zero-temperature critical line within the chirally symmetric phase, see fig. The free 
energy of the solution Zoo = 00 is 

Jl||(Zoo = 00) 



This result is exact. For the nontrivial solution we use 

b 



roo o 

/ du Vyp + i^ - -bn^ . (c.i) 

Jo 2 



^3/2 

as an expansion parameter in the ansatz for Zoo, 



(C.2) 



Zoo ^ SQ/ e + Qse^ • (C.3) 

Our result will show that along the critical line e is indeed small. Inserting this ansatz into 
( |4.6| ), expanding the right-hand side of this equation up to third order in e, and comparing 
the coefficients order by order yields 

Qi = ^^^^^, Q2 = ?Qf/^-|Qr^^^^^- (C.4) 

x/TT 4 2 x/TT 
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In order to find the transition between the two solutions we need to find the zero of AO, 



which is the difference in the corresponding free energies. Inserting eq. (C.3) into the free 
energy ([4.7|), subtracting (C.l) from the result and expanding up to second order in e yields 



AO , 3 _ 2 



with 



3^5/2 , r(^)r(f) 
2 



^1/2 r- 

Qi Vvr 



This yields at the critical line AO = 0, 

3 r~9 



4Q3 V l6Qi 7QsQl 



5/2 



0.0950977, 



(C.5) 



(C.6) 



(C.7) 



which is the result ( [4.10 ). Consequently, along the entire critical line e is constant and 
much smaller than 1, which validates our approximation a posteriori. 



D. Analytic approximation for zero-temperature chiral transition at small 
magnetic field 

Here we derive an analytic approximation for the T = chiral phase transition line at 
small values of the magnetic field b. This is not only a check for our numerical result but 
will also help to gain further insight into IMC. 

For the chirally broken phase we first solve eqs. (|3l0|) up to 0(6^), 



62 



uo ~ noo + 



cot ■ 



TT 

16 



V '^uoo ) 



3P2 
Pi 



(D.la) 



(D.lb) 



where the numbers Pi and P2 are defined in eqs. (|3.19 ) and where uqo = uo{b = 0) 



(2Pi)2/£2, The coefficient of the iJ.'^b'^ term in eq. ( p. lb ) is negative which implies that uq 



may in fact decrease as a function of b at fixed fi. Since uq is the location of the tip of the 
joined D8-branes and thus is proportional to the constituent quark mass, IMC manifests 
itself not only by a symmetry restoration but also more directly by a decreasing chiral 
condensate for certain chemical potentials. This decrease is only possible in the presence of 
a supercurrent, i.e., an anisotropic chiral condensate, which introduces the /i-dependence 
of Uq. (The dominant reason for IMC in the phase diagram is however the form of the 
free energy difference, as discussed in the main part of the paper, see also the following 
approximations.) Our numerical results show that for larger magnetic fields the value of 
Uq increases again for all fj, as expected from MC at weak coupling and converges to the 
value ( ^1T|) . 

Although we are interested in the free energy difference, it is instructive to consider 
the free energies of broken and symmetric phases separately. To obtain a finite result for 
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Figure 11: Comparison between tlie analytic approximation of the chiral phase transition (dashed 
Unes) with the numerical result (solid line) at zero temperature. The two dashed lines approximate 
the transition between the broken phase and the "higher Landau level" (hLL) phase and between 
the broken phase and the "lowest Landau level" (LLL) phase for small b. We have (unphysically) 
extended the phase transition line between the broken and the LLL phase into the parameter region 
below the (dotted) "Landau level" transition line. This is firstly a check for our approximation and 
secondly shows that IMC is strongest in the LLL phase, with the "higher Landau levels" working 
against it. 



the separate energies we subtract the vacuum contribution ^\\{t = fi = 0). Inserting eqs. 
( |D.1| ) into the free energy ( |3.14| ) yields, up to 0{b'^), 

-V 7 2 (2 16 

For the chirally symmetric phase we distinguish between the "lowest Landau level" (LLL) 
phase and the "higher Landau level" (hLL) phases (see sec. [4.1| ). The free energies of both 
phases are obtained from eq. (|4.7|). We find 



with Qi and Q3 defined in the previous appendix. To obtain Q^^^ we have used the 
expansion ( C.3 ). 

Consequently, the differences in free energies A^l = — flu between the broken phase 
and the two symmetric phases have the forms 

AOhLL ^ /i(^) - /2(m)6' , AQlll ^ const - ^fiH + ^(/i)^' , (D.4) 

where the constant and the functions /i(^), /2(a*)) /3(/^) can easily be read off from the 
previous equations. We plot the zeros of these free energy differences in the b-fi plane and 



compare them to the full numerical result in fig. 11 
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In the relevant regime of chemical potentials we have > for ^ < /ic(t = 6 = 0) 

and /2(/i) > for /i > 0.016/£^. This shows that if we start from the broken phase at 
6 = for 0.016/£^ ^ M < f^c{t = 6 = 0) and increase b at fixed n, chiral symmetry can 
be restored, provided that the sign change of AJI^ll occurs in a regime where our small-6 
expansion is valid. The comparison with the full result shows that our expansion is indeed 
a very good approximation to the full result, leading to IMC in the hLL regime. 

The transition between the chirally broken phase and the LLL phase occurs at relatively 
large magnetic fields where our expansion is a less accurate, however qualitatively still 
reliable, approximation. In the figure we have extended the phase transition line between 
the broken and the LLL phase into the region where the hLL phase is the ground state. 
This extended line is in very good agreement with our approximation and shows that IMC 
is strongly dominated by the LLL with the hLLs working against it. 
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